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Abstract 

We show that the low-momentum behavior of the Landau-gauge gluon and ghost propagators changes on the lattice if the gauge- 
fixing procedure favors Gribov copies with an exceptionally small lowest non-trivial eigenvalue of the Faddeev-Popov (FP) operator. 
Compared to random copies, the ghost propagator below 1 GeV grows stronger towards zero momentum on Gribov copies with 
^vq a small lowest-lying FP eigenvalue, whereas the gluon propagator is more suppressed below 0.2 GeV. Above these momenta no 
i—\ dependence on Gribov copies is seen. Qualitatively, our data thus resembles the change of the gluon and ghost dressing functions, 
C> Z and J, with the boundary condition on 7(0), put forward by Fischer, Maas and Pawlowski [Annals Phys. 324 (2009) 2408]. 
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1. Introduction 

Lattice calculations of the Landau-gauge gluon, ghost and 
quark propagators have attracted quite some interest during the 
last 15 years. Staunch supporters of pure lattice QCD (LQCD) 
may wonder about the enthusiasm with which such calcula- 
tions have been performed and discussed in the past, in par- 
ticular, as LQCD comes with the distinct advantage that one 
does not need to fix a gauge. This holds true, however, only 
as long as one is interested in gauge-invariant quantities. But 
besides LQCD there are also other (sometimes better suited) 
frameworks to tackle nonperturbative problems of QCD, and 
these require the exact knowledge of QCD's elementary two 
and three-point functions in Landau or other gauges. 

Two continuum functional methods one has to mention here 
are the efforts to solve the infinite tower of Dyson-Schwinger 
equations (DSEs) of QCD or, likewise, the corresponding Func- 
tional Renormalization Group Equations (FRGEs) (see, e.g., 
the reviews |[l}|6) and references therein). Both these meth- 
ods imply fixing a gauge (and often the Landau gauge is chosen 
for simplicity), but more importantly, these methods also re- 
quire a truncation of the infinite system of equations to enable 
finding a numerical solution. These truncations are a potential 
source of error, which why corresponding (volume and contin- 
uum extrapolated) lattice results are so essential to render these 
truncations harmless or to even substitute parts of the DSE (or 
FRGE) solutions by (interpolated) nonperturbative data. 

In what concerns the Landau-gauge gluon and ghost propa- 
gators, lattice results have helped much to improve truncations 
over the years. Currently, the continuum and lattice results 
overlap for a wide range of momenta, showing nice consis- 
tency among the so different approaches to QCD. Admittedly, 



the currently used truncations are still not perfect, as seen, for 
example, for the gluon propagator whose DSE solutions dif- 
fer from the corresponding lattice or FRGE results in the in- 
termediate momentum regime (i.e., for momenta 0.5 - 3 GeV), 
whereas FRGE and lattice results agree much better there (see, 
e.g., Fig. 2 in JT]). But this situation will certainly improve, as 
it did in the past. 

Another regime that remains to be fully settled yet is the low 
(infrared) momentum regime. About the infrared behavior of 
the gluon and ghost propagators in Landau gauge there has been 
much dissent in the community and it is difficult to assess on 
the lattice also. Currently, all lattice studies agree upon a gluon 
propagator and ghost dressing function which are (most likely) 
finite in the zero-momentum limit (see, e.g., EHT6|^] DSE 
and FRGE studies fl7l - l20ll . on the other hand, assert that this 
infrared behavior is not unique, but depends on an additional 
(boundary) condition on the ghost dressing function at zero mo- 
mentum, 7(0). Explicitly, in Refs. {19\ 1201 . it is shown that 
for / (0) = one finds the so-called scaling behavior for the 
gluon and ghost propagators at low momentum, as it was first 
found in [21 ], while for finite J(0), one finds a family of decou- 
pling solutions for the DSEs and FRGEs, in qualitative agree- 
ment with DSE solutions proposed in the studies of Refs. 11221 - 
l27ll . and with lattice results. For momenta above 1 GeV both 
types of solutions are practically indistinguishable. We inter- 
pret this ambiguity in the infrared as a remnant of the Gribov 
ambiguity of the Landau gauge condition, which is lifted by 
fixing /"'(O) to a constant. 

In this letter we will show that this one-parameter family of 
decoupling solutions can also be seen on the lattice, at least to 
the extent computational resources allow. After outlining some 
technical details in the next section and a discussion about the 
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distribution of the lowest non-trivial eigenvalue of the Faddeev- 
Popov (FP) operator, A\, on different Gribov copies (Sect.[3]l, 
we will demonstrate in Sect.|4]that the decoupling-like behav- 
ior of the lattice gluon and ghost propagators can be changed by 
a (yet simple-minded) minimization of A\. Changes take place 
only in the low-momentum regime, but interestingly in a simi- 
lar manner as shown in [19], where the condition on 7(0) was 
changed. Specifically, we show that on Gribov copies with very 
small A\ the ghost dressing function at low momenta is more 
enhanced than on random copies, while the gluon propagator is 
more suppressed. 

We still find Gribov copies by a maximization of the lat- 
tice Landau-gauge functional, but we are not interested in find- 
ing Gribov copies with a large (local) gauge-functional value, 
but on copies with very small A\, irrespective of the functional 
value. Only then we see a stronger rising ghost dressing func- 
tion and a more suppressed gluon propagator at low momen- 
tum. On Gribov copies with large gauge-functional values we 
see both being suppressed, consistent with what was found in 
the past Il8l [28H3TT1 . 

We should mention here that similar effects were seen for the 
B-gauges by Maas [32 1. For these gauges, one selects Gribov 
copies based on the ratio of the ghost dressing function at a 
small and a large lattice momentum on a particular copy. By 
construction the ghost dressing function in these gauges is then 
clearly enhanced or suppressed at low momenta. It remains to 
be seen if corresponding effects become clear also for the gluon 
propagator. 



2. Simulation details 

Our study is based on 60 thermalized gauge field configu- 
rations, generated with the usual heatbath thermalization and 
Wilson's plaquette action for SU(2) lattice gauge theory. The 
lattice size is 56 4 and the coupling parameter /3 = 2.3. To 
reduce autocorrelations, configurations are separated by 2000 
thermalization steps, each involving four over-relaxation and 
one heatbath step. For every configuration we generate at least 
^Vcopy = 210 gauge-fixed (Gribov) copies, all fixed to lattice 
Landau gauge using an optimally-tuned over-relaxation algo- 
rithm for the gauge fixing. To ensure these Gribov copies are all 
distinct, the gauge-fixing algorithm always starts from a random 
gauge transformation field. Interestingly, for all these 60 x 210 
gauge-fixing attempts we only found one Gribov copy twice. 

With the help of PARPACK ED we then calculate the low- 
est three (non-trivial) eigenvalues < A\ < A% < /I3 of the 
Faddeev-Popov (FP) operator of every Gribov copy, and clas- 
sify them by their value for A 1 . The Gribov copy with lowest Ai 
(considered for each configuration separately) is labeled lowest 
copy (£c). The first generated copy, irrespective of A\, gets the 
label first copy (fc). It represents an arbitrary (random) Gribov 
copy of a configuration. To compare with former lattice stud- 
ies on the problem of Gribov copies we also reintroduce the 
label best copy (be). It refers to that copy with the best (largest) 



gauge functional value 

Fu[ ^ = ^ZZ^Tr^C/^. (1) 

X [1=1 

for a particular gauge configuration U = {U x f,}. Here g = {g x } 
denotes one of the many gauge transformation fields fixing U 
to Landau gauge. 

On those three sets of Gribov copies we calculate the SU(2) 
gluon and ghost propagators following standard recipes. That 
is, the gluon propagator is calculated for every lattice momen- 
tum using a fast Fourier transformation and the ghost propaga- 
tor by using the plane-wave method for selected momenta. To 
accelerate the latter we use the preconditioned conjugate gra- 
dient algorithm of (8). As a by-product of this calculation we 
also obtain the renormalization constant, Z\, of the ghost-ghost- 
gluon (gh-gl) vertex in Landau gauge for zero incoming gluon 
momentum. For more details on lattice Landau gauge and the 
calculation of the propagators and Z\ the reader may refer to 
Refs. [8 34-36 1 and references therein. 

When quoting momenta in physical units we adopt the usual 
definition ap^k^) = 2 sin (nk^ / L^) [with k^ e (-Z^/2, L M /2] 
and Lft = 56], assume for the string tension -\fo- = 440 MeV and 
use era 2 = 0.145 for /3 = 2.3 from Ref. (37), where a denotes 
the lattice spacing. 

3. Distribution of ii 

Before comparing the propagator data for the different types 
of Gribov copies, it is instructive to look at the distribution of 
A\ on all copies first. In Fig. [T] we show this eigenvalue distri- 
bution (in lattice units) for N cp = 210 Gribov copies. There, the 
big panel shows it separately for each of the 60 gauge configura- 
tions and the small panel (on top) for all configurations together 
as a histogram. One sees that for most of the copies Ai takes val- 
ues between 0.5 x 1 0~ 3 and 1 .9 x 1 0~ 3 , mostly between 1 .5 x 1 0~ 3 
and 1 .7 x 10~ 3 , but for some configurations there are also copies 
with an exceptionally small value for A\ , a value [Ai < 10~ 4 ) far 
below the values found for the other copies. With our simple 
(brute-force) approach we are rather limited in finding more of 
these exceptional copies. The gauge-fixing and calculation of 
eigenvalues on a 56 4 lattice is computational quite demanding, 
and a more sophisticated gauge-fixing algorithm — one which 
would automatically select that Gribov copy with the smallest 
(or at least with small) A\ — does not exist. But it would be in- 
teresting to know if for each configuration a Gribov copy with 
such an exceptionally small A\ exists. 

For a few configurations we generated more than 210 Gribov 
copies. These allow us now to have a closer look at the dis- 
tribution of Ay and to demonstrate that this number of copies 
is sufficient to resemble the distribution's shape for our lattice 
parameters (J3 = 2.3, 56 4 ). 

Typical snapshots of this distribution for different N cp are 
shown in Fig. [2] from left to right for three arbitrarily selected 
configurations. There, the middle panels illustrate how the 
range of Ai -values gets populated when increasing N cp ; shown 
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Figure 1: Distribution of A\ (in lattice units) for 210 Gribov copies, shown 
separately for 60 thermalized gauge configurations (ordinate). Full red (white) 
circles mark the eigenvalue on tc [fc) copies. The small top panel shows as a 
histogram for A\ including all values found. 



Gribov copies (N cp = 420 or 500) for the respective config- 
uration. In former lattice studies on the infrared behavior of 
the gluon and ghost propagators (e.g., lfT2"l [131 [3 , 38 1), gauge- 
fixing algorithms were often designed to find copies with com- 
parably large Fy[g], in the hope these copies are closer to the 
fundamental modular region than a set of random copies can 
be. Gribov copies which globally maximize the gauge func- 
tional are elements of the fundamental modular region, and in 
the continuum the common boundary of this region and the Gri- 
bov horizon (the set of Gribov copies with A\ = 0) is expected 
to dominate the path integral in the thermodynamic limit [39]- 
HTI . One might thus expect that a small value for A t is corre- 
lated to a large value for Fu[g]. Looking at the distributions 
for A\ and F u[g] (Fig. [2]), we find, however, there is no obvi- 
ous correlation between them. There are copies with small A\ 
and small Fu[g], but at the same time there are also copies with 
small A\ and large F v [g], and vice versa. 



4. Ai and the gluon and ghost propagator 

Next we look at the correlation of A\ and the gluon and ghost 
propagator at low momenta. For the gluon propagator D(k) [or 
its dressing function Z(k) - p 2 D(k)] there is actually no direct 
relationship between its momentum dependence and the spec- 
trum of the FP operator, besides that it is gauge-dependent and 
so may change if the lattice Landau gauge condition is supple- 
mented by a condition on A\ . 

The ghost propagator G(k) [or its dressing function J(k) = 
p 2 G(k)], on the other hand, should be affected stronger. Given 
its spectral decomposition in SU(2), 
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are distributions for jV cp = 70, 210, 420 and 500 Gribov copies. 
The top panels show the corresponding histograms, filled with 
the respective symbol color used in the middle panels. To ease 
the comparison, these histograms are all normalized with re- 
spect to jV C p = 420. 

From these histograms we see that the individual 
(configuration-wise) distributions of Ai are asymmetric 
and negatively skewed, similar to what we have just seen for 
the overall distribution in Fig.[T| Moreover, we see that at least 
200 Gribov copies are needed to reach at an approximate shape 
for the distribution. This number of copies seems to be also 
sufficient (for the given lattice parameters) to find a copy with 
a very small A\, even though it is unlikely that this copy is the 
one with the smallest A\ overall. 

We are thus in a good position to analyze the relation of a 
very-low A\ and the low-momentum behavior of the gluon and 
ghost propagator. 

The attentive reader may have noticed the bottom panels of 
Fig. [2] These show correlation plots of A\ versus the gauge 
functional Fy[g], always for the largest available number of 



it is plausible that Gribov copies with very small Ai yield larger 
values for G than copies where Ai is comparably large, possi- 
bly even configuration- wise. Though this is not as simple as 
it seems at first sight, because A t and the corresponding eigen- 
function <t>i(&:) contribute only a minor fraction to G(k) and this 
fraction even shrinks the larger the lattice momentum a 2 p 2 (k) 
(see Ref. [42 1, in particular Fig. 6 therein). From [43 ] (Fig. 1) 
we learn that this fraction seems to increase with volume. 

For a few gauge configurations we have data for the gluon 
and ghost propagator for all Gribov copies. It allows us to 
check if there is a direct correlation between A\ and the gluon 
and ghost propagator at low momenta, comparing different Gri- 
bov copies of the same configuration. Looking at this data re- 
veals, however, no immediate correlation: A Gribov copy with 
a smaller A\ not necessarily yields a larger ghost propagator 
than a copy with a somewhat larger A \ . Similar we find for the 
gluon propagator. 

Nonetheless, when looking at these correlations more 
broadly, that is on average for the whole gauge ensemble, we 
see a clear trend: Gribov copies with very small A\ tend to yield 
a larger ghost propagator at low momentum, than copies with 
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Figure 2: Top panels: histogram of the eigenvalue distribution of A[ for different numbers /V cp of Gribov copies. From left to right, each panel is for one gauge 
configuration. All histograms are normalized such that the sum of bars is one for Nop = 420 copies. The middle panels show the corresponding scattering of A[ for 
the different N cp , and the lower panels show scatter plots of Ai versus the gauge functional values Fu[g\ for the respective largest N cp . 



large A\. This anti-correlation is smaller the larger the momen- 
tum, and is also considerably smaller for the gluon propagator 
at same momentum. 

We clearly see both these correlations when comparing aver- 
aged lattice data as shown in Fig. [3] There the left panels com- 
pare the momentum dependence of the SU(2) ghost dressing 
function (top) and the gluon propagator (bottom) as obtained 
on first, lowest and best copies of our gauge ensemble (refer 
to Sect.|2]for this classification). The right panels confront the 
corresponding data for the strong coupling constant, 
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here in the Minimal MOM scheme for SU(2) Landau gauge 
B41 . and of the inverse renormalization constant Z~' of the 
gh-gl vertex with zero incoming gluon momentum. 

Looking first at the first-copy (fc) data points (black open 
circles in Fig. [3]), we see these behave as expected IT51 [TT1 fOl 
|30j[3T|]: The ghost dressing function and the gluon propaga- 
tor increase with decreasing momentum and tend to reach a 
plateau at very low momentum, while ar™ 2 ) grows with de- 
creasing momentum down to about a 2 p 2 = 0.2 (p « 0.5 GeV) 
below which it falls off again with momentum. Also Z~ l is 
found as expected: Zj" 1 — > 1 for large momenta, as it should, 
and for intermediate momenta we see its characteristic hump 
HUES] 02. 

But the most interesting data in Fig.[3]is the lowest-copy (tc) 
data (red full squares). The ic data points, for example, for the 
ghost dressing function clearly deviate upwards from the re- 
spective fc data below a 2 p 2 = 0.5, while for the gluon propaga- 
tor these systematically deviate downwards below a 2 p 2 - 0.03. 
These deviations then yield a coupling, a^(2y whose running 
is slightly upwards shifted below a 2 p 2 = 0.5. Admittedly, com- 
pared to the ghost dressing function, the effect for the gluon 



propagator is small, but we find that it becomes more and more 
pronounced when increasing statistics. Looking, for instance, 
at half the number of configurations the effect is less signifi- 
cant than for the full set. So we think the effects we find for 
both propagators and a^ (2) are not statistical artifacts, but re- 
sult from our selection of Gribov copies with very small A\ . 

For completeness, we also show the corresponding best-copy 
(be) data in Fig. [3] (blue crosses). Comparing this with the 
fc data, we find the expected suppression of the be data for 
the ghost and gluon propagator at low momentum which is then 
also seen for ^supy Similar was seen for the be data in lfT2ll30l . 
although there the effect was even bigger when using the FSA 
gauge-fixing, a special combination of the simulated annealing 
and over-relaxation algorithm for gauge-fixing and Z(2) flips. 

For all three types of Gribov copies, we find that Zj" 1 remains 
almost unaffected (see lower right panel of Fig.|3]l. There are 
small upward shifts for the Ic data at small and at large momen- 
tum (there also for the be data), but these shifts are all within 
(statistical) errors, besides at lowest momentum. This certainly 
deserves further study, because at large momenta Z\ — > 1 for all 
types of copies [44. 46 1, while at small momenta Z\ is expected 
to approach one or at least a value close to one fl7l . 



5. Summary and conclusion 

We have demonstrated that the low-momentum behavior of 
the Landau-gauge gluon and ghost propagators can be changed 
on the lattice by an additional constraint on the lowest-lying 
(non-trivial) eigenvalue, A\, of the FP operator. If the lattice 
Landau gauge fixing is tuned to find Gribov copies with a very 
small A\, the ghost dressing function gets more enhanced to- 
wards the infrared momentum limit and the gluon propagator 
more suppressed. The combined effect then yields a coupling 
constant ff™ 2 ) whose running is slightly upwards shifted for 




Figure 3: Ghost dressing function (top left), gluon propagator (bottom left), the coupling aSjjj™ ( t0 P r 'S nt ) ant ^ me inverse renormalization constant of the gh-gl- 
vertex (bottom right); all versus lattice momentum squared, in lattice units and for SU(2), and no renormalization has been applied. Open circles (filled squares, 
crosses) refer to data on fc (fc, be) copies. The lower left panel also shows a zoomed-in plot to improve visibility of the low-momentum region. Note the different 
scales for the x-axis (log and linear) in the left panels and the corresponding physical momenta above the top panels. 



momenta below 0.7 GeV. In contrast, Gribov copies with a large 
gauge functional value tend to suppress instead both propaga- 
tors and ar™ 2 ) at l° w momentum. 

These Gribov-copy effects are not noticeable at momenta 
above 1 GeV. Also, these effects cause only a quantitative but 
no qualitative change of both these propagators at low momen- 
tum. Their modified momentum dependences are still of de- 
coupling type. But interestingly, the change we find with A\ 
looks very much alike the change of the gluon and ghost dress- 
ing functions, Z and J, with the (boundary) condition on /(0) as 
found in [ 19 Looking at Fig. 5 therein, one sees that if 7(0) is 
fixed to a larger value the ghost dressing function below 1 GeV 
is enhanced, while the gluon dressing function below 0.2 GeV 
is decreased. One further sees that the effect is stronger for the 
ghost than for the gluon dressing function (compared at same 



2 Note that in 1 19 | the ghost dressing function is denoted G. 



momentum), and that the deviation of the strong coupling con- 
stant appears below and slightly above the position of its maxi- 
mum at a momentum of about 0.5 GeV. 

This is exactly what we find from our Ic data, suggesting that 
the constraint on J(0), used in [19| to lift the Gribov ambigu- 
ity of Landau gauge, can be mimicked on the lattice by sup- 
plementing the lattice Landau-gauge condition by a constraint 
on Ai, 

Strictly speaking, our results are for SU(2) only, but similar 
should be expected for SU(3). The propagators for these gauge 
groups differ only little IfTTl l47ll . SU(2) and the rather coarse 
56 4 lattice (J3 = 2.3) was chosen here on purpose, because it has 
allowed us to get (with reasonable computer time) data for the 
gluon propagator and the ghost dressing function where their 
momentum dependence starts to become flat, and this for a suf- 
ficient number of Gribov copies (at least 200 per config.) such 
that a correlation between A\ and both propagators can be seen. 
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For a future study, we suggest to cross-check our results at 
lower momenta, in particular, for the gluon propagator whose 
low-momentum behavior changes only little, if one constrains 
X\, or 7(0) in the continuum. Also 3-point functions of gluon 
and ghost fields should be checked, and what changes when 
adding fermions. Note that larger lattice sizes and smaller lat- 
tice spacings will perhaps also affect the low-momentum be- 
havior on any finite lattice [16, 36 1 . Such calculations will thus 
become expensive very quickly — especially, if Gribov copies 
with exceptionally small A\ still have to be found by chance. 
Algorithmic improvements would therefore be quite helpful. 
For example, a gauge-fixing algorithm that automatically se- 
lects Gribov copies with very small A\ . 

Let us comment on the "rigorous bounds" introduced in 
Il43l |48 1 to control the infinite-volume extrapolations of lattice 
data for the gluon and ghost propagator. Since these bounds 
are composed of the lattice gluon field (for the gluon propaga- 
tor), and of A\ and the eigenmode <5i (for the ghost propagator), 
these will suffer similar Gribov problems as seen for the propa- 
gators here. For a particular lattice implementation of Landau- 
gauge Yang-Mills theory and a particular selection of Gribov 
copies these bounds will certainly constrain the lattice data, but 
one should keep in mind these bounds may change when using 
others. 

With the current lattice implementation of Landau-gauge 
Yang-Mills theory, one seems to be able to realize different 
decoupling-type solutions for the propagators by supplement- 
ing the lattice Landau-gauge condition by an additional con- 
straint, as, e.g., the one here or in (32), but it will perhaps be 
impossible to find the scaling solution on a finite lattice. Even 
if we were able to find for every gauge configuration that Gri- 
bov copy with the smallest X\, this eigenvalue and all the others 
would still be non-zero, i.e., this copy is not on the Gribov hori- 
zon. It remains to be seen if for such an approach the scaling 
behavior can then appear in the limit of infinite volume and zero 
lattice spacing. 
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